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In the recent years many researches were performed about graphene, as we can see frorrP®^. 
Graphene is always considered a half metal or a zero gap semiconductor. In the last year new 
experiments were done about graphene on boron nitrid^, and in particular in a work they obtained 
an insulating behaviour and a power law for conductivity at the charge neutrality point (CNP) or 
Dirac point (DP)P. In this work we will explore the electron-hole asymmetry and the disorder in 
graphene. In particular we analyze the asymmetry in the ordered phase, considering thermal doping 
at half filling and with a non zero chemical potential. After that we use the self consistent Born 
approximation for calculating the self energy correction from impurity scattering. 

Furthermore we shall analyze the electron-phonon interaction, using the deformation potential model 
for obtaining the temperature dependence for the compressibility. Finally we will speculate about 
temperature dependence of conductivity. 


I. INTRODUCTION 

Graphene is a mono layer material with a honeycomb 
lattice made by carbon atoms and it is a relative of 
graphite. It can be considered the first two dimensional 
material, since its discovery in 2004^. In fact Geim and 
Novoselov discovered how to stabilize a two dimensional 
material by means of a substrate, against the general 
opinion motivated by the Mermin-Wagner theoreirP. In 
the last years many scientific works were dedicated to 
Graphene likj ^ wh ere we find many of the electronic 
properties ancP^^^ about experimental properties, be¬ 
cause graphene is a zero gap semiconductor or a half 
metal with a very peculiar band structure near the Fermi 
leveP^ with a linear dispersion instead of a parabolic one 
like many other semiconductors. This properties bring 
with them many peculiar and fascinating characteristics 
that can be used in applicationJi]. Up to 2013 All the 
studies were performed on graphene on Si/Si 02 but it 
was pointed that this kind of substrate can modify the 
properties of graphene with a modification of the Fermi 
velocity and the band structur^i^. Just in the last two 
years new experiments were performed on a substrate of 
boron nitride, a two dimensional material with atoms of 
boron and nitride arranged in a honeycomb structure but 
with a semiconducting behaviour Furthermore all previ¬ 
ous experiments where performed with a single gate set 
up, checking just the doping of the system. In a recent 
experimeniP They performed measures on graphene de¬ 
posed on boron nitride with a double gate experiments 
which allows them variations of the number of particles 
and also the chemical potential. Furthermore new tech¬ 
nology advances brought to the market new organic gates 
with an increase efficiency. They discovered an insulat¬ 
ing behaviour near the CNP and a very strong temper¬ 
ature dependence of the conductivity as a power law. 
In this work we aim to derive some important proper¬ 
ties of graphene near the Dirac point and the effects of 
the electron-hole asymmetry. In particular in the sec¬ 
ond section we shall analyze the e-h asymmetry on the 


density of states and electronic properties. In the third 
section we will see the thermal doping in the presence 
of the asymmetry in the non interacting system. In the 
fourth section we will see the role of the static and dy¬ 
namical disorder by means of the Self Consistent Born 
Approximation (SCBA) for impurity scattering and the 
electron-phonon interaction and the derivation of the self 
energy and the interacting DOS with the spectral proper¬ 
ties. Moreover we will analyze the effects of the disorder 
in the Dirac point and finally in the last section we shall 
see the compressibility of the system that may be linked 
to a quantitative explanation of the anomal conductivity. 

II. LOW ENERGY ELECTRONIC MODEL AND 
E-H ASYMMETRY 

In this section we are going to introduce our low energy 
model for describing the electron-hole asymmetry. As we 
have seen so far the graphene’s band structure near the 
Fermi level presents itself like two cones touching in one 
point and the density of states is a linear function of 
the energy with a different sign for electrons and holes 
as we can see frorrP^. In particular we will calculate a 
correction to the bare density of state (DOS) in the Dirac 
point (DP), obtaining a non linear behavior. In fact if 
we consider a Taylor expansion of the band structure up 
to the second order, as obtained bjlH, we have: 

E±{q,9) =3t ±Vfq- A{e)q‘^ (I) 

where A{9) = ± cos[30]. If we use the Jacobian 

of the transformation and the lower and upper values 
of q and 9 the DOS is p{E) = d9 6{E - 

E±{q, 9))qdq, where = |ai x 02 ! is the volume of the 
unitary cell. The calculation can be simplified by means 
of a shift of the energy variable as e = E — 3t =0 and 
dividing the integral between electrons and holes, where 
we use a plus sign at the right of the DP and a minus one 
at the left. We notice that in the limit oft —)■ 0 we obtain 
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again the expression of a symmetric DOS. Furthermore 
we shall write the DOS in a more compact form as: 

p±(e) = (a|e|+sgn(e)/3e^) (2) 

where a = and P = Using a 

value of t = 0.2t as reported in literature' , we obtain 
the numerical values (in hopping units) a ~ 0.368 and 
P ~ 0.221. 


III. THERMAL DOPING IN 
NON-INTERACTING SYSTEM 

In this section we are going to derive the thermal dop¬ 
ing in the non interacting graphene. As we have seen so 
far the new experiments concerning the double gate dop¬ 
ing in nano-structures are suitable to control not only the 
number of particles but also the chemical potential, so it 
might be important to have an analytical expression for 
that. As previously pointed out by Guinea et al.^ we 
have thermal doping just with electron-hole asymmetry. 
In particular we may have three different cases, distin¬ 
guished by the values of the chemical potential: 

1. The DP with fj, = 0, analyzing the number of par¬ 
ticles as a function of temperature; 

2. The case with a non zero chemical potential /i ^ 0 
which can be split in two sub-cases of low and high 
temperature: 

• the low temperature limit \fj,\ < T that is equal 
to the DP thermal doping plus some correc¬ 
tions derived by an expansion of the Fermi- 
Dirac function; 

• the high temperature limit \fi\ > T where it is 
possible using the Sommerfeld expansion with 
properly used cut-off; 


A. Thermal doping at the Dirac point 

If we fix the chemical potential in a double gate experi¬ 
ment we are able to use the symmetry of the Fermi-Dirac 
function for evaluating the thermal doping. In fact divid¬ 
ing the DOS in a symmetric part (first order in e) and a 
non symmetric part (second order in e) adding plus and 
minus ^ to the expression for the number of particles we 
have: 


^ ^ + oo / \ 

n{^i,T) = - j {ps{e)+ 

(3) 

where Ps{f) = Oi\e\ e Pas{f) = sgn(e)/3e^. We notice that 
the first integral is zero, because the function is odd in 
an even interval, so we shall solve just the second one. 
Using a simple substitution ^ = x, we obtain the final 


expression with a power law function for the temperature 
in the DP: 

n{p = 0, T) - n{p = 0, T = 0) = 8 T^pi 2 (4) 

with I 2 = x^(l — tanhxjdx = |C(3), with the nu¬ 
merical value: 

n(T) -no = 3.61PT^ (5) 

where the number of particles for a half filling not doped 
system is rig = 2 with a maximum equal to 4. We notice 
that for /3 —1 0, i.e. t —0 we do not have a temperature 
dependence of the number of particles, because it is a 
direct consequence of the electron-hole asymmetry of the 
density of states. 


B. Non-zero chemical potential 


The low temperature regime is characterized hy p ^ 0 
and p ^ T and with these conditions we can use the 
Sommerfeld expa nsion. Following the technique reported 
in literatur d^^bs i -^i-jte the number of particles 

with both contributions from electrons and holes as: 

/ O /•+00 

p-{e)f{e)de+ / p+{e)f{e)de (6) 

-00 J 0 


where p- and p+ refer to the density of states at left 
and right of the DP, respectively. In the limit of low 
temperature and expanding the integrals in a series: 


n{p,T) -no = sgn{p) 
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T‘^{a + 2Pp) 


(7) 


where we see that is a small corrections to the thermal 
doping induced by /3 and for T —>■ 0 we obtain the integral 
of the density of states of the system. The final expression 
in the high temperature limit is: 


n{p,T) - no = 

3MPT^ + p (2aT/i -F 2pT'^l2) - p^ {alg + 2pTh) 

( 8 ) 


where we notice a symmetry behaviour for /i —>■ 0, with 
an expression given again by eq.([^. We see the eq.Q 
as a function of temperature and the chemical potential 
in fig.([^ with a comparison with the exact numerical 
result. The arrow set the limit of validity of the analytical 
solution. 


IV. ROLE OF STATIG AND DYNAMIGAL 
DISORDER 

A. SGBA -|- el-ph interaction 

In this section we are going to see the effects of the in¬ 
teractions on the system. In fact we have seen all the 
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FIG. 1. Particle number as a function of the chemical po¬ 
tential (fixed temperature) and Temperature (Fixed chemical 
potential). In both pictures kb = 1, the energy scale is in 
hopping units t = 0.2t and everything is shifted by no = 2 


previous results with a zero temperature and without 
interactions. In particular following previous works we 
can consider the effects of impurities by mean s of the 
Self consistent Born approximation (SCBA^P^Hland also 
the combination of the SCBA and the electron-phonon 
interactiorP^. In particular if we consider a potential 
for short-range scatterers like in AndcJ^^jWith the Dyson 
equation and the SCBA the Green’s function will be a 
diagonal matrix, indipendent from k and the self energy 
can be written as a product of the averaged potential 
with the averaged Green function with a Gaussian disor¬ 
der, so: 


E(e) 

ma^ 




de 


P+(e') 

e- e' - E(e) 


(9) 


where the bosonic propagator is: 

D{q,iUm) = - 




-h 


and the Fermionic propagator: 
G(k-|-q,ia;„) = 


I 


iuJn + M - ek - s(fa;„, k) 


( 11 ) 


( 12 ) 


and and Vm are respectively the Fermionic and 
Bosonic frequencies as reported in MahaiP^. 

Considering the approximation of classical phonons 
Tiujq ^ kbT we can still write the self energy as indepen¬ 
dent from k, like in the SCBA from impurity scattering, 
with a propagator given by: 

2 

D(q, c::i djji Q (m) 

UJq 

and so the Self energy, a totally isotropic function in the 
IBz, is 

S(ia;„,k)=ry^ [ duj2- --|-q, iw„) 

q Jo w 

(14) 

where Q!^F(q, w) = /|q|(5(a; — Wq) is the Eliashberg func¬ 
tion and / is linked to the deformation potential as re¬ 
ported by MaharPSl. In this way we obtain: 

S(za;„,k) ~ cr^G^Q(,(iw„) (15) 

where the variance of the disorder is the same of the 
variance obtained for the impurity scattering. If we write 
the variance by means of the microscopic parameters of 
the e-ph interaction, with the deformation potential we 
have: 


2TI 

- = T - 

^Vs yPmVl 


(16) 


with D deformation constant, Vs the sound velocity, pm 
the density and V the volume of the unitary cell. Using 
the values reported in literatur^^ we obtain: 


t2 


I.3( 


kbT 

t 


) 


(17) 


which is the variance of the disorder for t = 0.2t and it 
is a linear function of the temperature. 


B. Energy dependent properties 


where A_|_ e A_ are integration cut-off, different because 
the DOS is asymmetric, cr^ ig the variance of the Gaussian 
disorder and rii the density of scatterers. Considering 
also the electron-phonon interaction we can write the self 
energy of the electron as: 


S(iU!n + = 

= ^ \gq\'^D{q, jzz„)G(k -|- q, iuJn + Wm.) 

^ m,q 


We have seen that we can obtain an expression of the 
self energy of the system by means of eq. (|^ . and finally 
we arrive to an expression for the self energy which in¬ 
cludes the effect of the e-h asymmetry: 


= a f A_ — A_]_ -b z In 


S(e) 

nicr2 

/ , , A2 -bA2 

+ /3 zG{z) - -± -= 


(2;-b A_)(z-A+) 


(18) 


( 10 ) 
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where we notice that using the substitution z = e — S(e) 
we obtain the self consistency and with the limit /3 —>■ 0, 
also = A_ and we obtain a similar expression like in 
AndcP^. 

If we neglect the self consistency in eq.( 181 we can obtain 
an analytical expression for the self energy of the system 
and analyzing the spectral properties. In particular using 
eq.([7l|) we can derive the expressions of the real (A) and 
imaginary (F) parts of the self energy: 


r(a;) = —niCT^7rsgn(a;)(atLi + 

A(w) = nia'^[{a + /3w)(A_ — A+ + 

''(v(w + A_)(a;-A+)j^ ^ 2 

that are expressions valid inside the energy between the 
cut-off. We notice the minus sign of the imaginary part 
of the self energy, needed because it is directly related 
to the density of states and again we obtain the sym¬ 
metric version imposing the symmetry as /3 —> 0. If we 
change the sign of the imaginary part we can compute 
the spectral DOS as: 


C. Dirac point properties 

We shall be able to analyze the disorder corrections to 
the Dirac point starting from eq.(18) at fixed energy,i.e. 


the DP as A = ReE = ui and so F = ImE. In this 
condition the real part is a correction to the electronic 
energies, a band renormalization, and the imaginary part 
is a correction to the density of states of the system. In 
this condition the self energy will be E = — zF and so 
z = zF. Remembering the expressions of the functions in 
eq. (691 we obtain: 


G(zF) = (A_ - A+)+ 
-F2 

zF In 


(zF-hA_)(zF-A+)_ 


(23) 


A2 4- 

F(iF) = iFG(zF)- ~ ^ + 


Now at fixed energy we can neglect the contribution of F 
inside the logarithm, because A+, A_ F, so: 


G(zF) ~ (A_ - A+)-hzFln 


F2 

aIa 


- I -. 


(24) 


1 

DOS{uj) = -[ / dep-{e)lm 

Ja. 


' — e — E(a;) 


( 20 ) 


dep+(e)Im- 


u! — e — E(a;) 


and so we obtain an expression for the interacting density 
of states: 


DOS{oj) = 

- ilm[(„(A_ - A+ + . log ^ 


)+ 


-I- /3z(A_ - A+ -I- z log 


(z-f A_)(z-A+) 


)- 


( 21 ) 




Ai -H A?., 


where z = w — A(a;) -I- zF(a;) with A(a;) = ReE(a;) and 
F(a;) = ImE(a;). Considering the different Riemann sur¬ 
faces we have the final expression: 


DOS{u) = 

-{(aF(w) -I- 2,5F(tLi)(a; — A(a;))) log + 

TT \/P2P3 

+ sgn(w)(20i - 02 - O 3 ) 

[a{oj - A) + ^{{uj - A{uj)f - F(a;)2)] + 

+ F(c.)^(A_-A+)} 


where the pi and 9i are functions of to as reported in the 
appendix. In the next section we are going to analyze 
what are the properties spectral properties in the DP. 


and the final equations will be: 

= [(a -I- Z/3F) ( A_ — A_|_ -1- zF In 




A_A+ 


-/3 


A^ -t A2 


+ 1 


(25) 


which is an explicit solution of eq.(18) at fixed energy 
w = A = ReE. 

Considering separately the contribution of the imaginary 
part, from eq.(25) we obtain: 


-F 

rigCr^ 


= F ( j3{A- — A+) -I- a In 


A.Aa 


(26) 


and in the same limit of A_, A-(_ ^ 1 we have an expo¬ 
nential small correction to the density of states which is 
different from zero at the DP: 


F = A/A+A_e"5)j(;(j+/5(A--A+)) 


(27) 


we notice that in the limit of /3 —>■ 0 the eq.(27) is the 
same of AndcP^, because A_ = A+. With a perfect anal¬ 
ogy we can define a renormalized disorder: 


l + /3a2(A_ - A+) 


(28) 


with a variance of the effective disorder that can be in¬ 
creased if j3 has a negative sign and it can be decreased 
if P has a positive sign, because A_ > A+ as we will see 
later. This is a totally new results and it is important its 
connection to both the static and dynamic disorder. In 
fact if we consider a system with impurities, combining 
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the effect with the e-h asymmetry we can increase the 
effect of the renormalization on the DOS and in a system 
with dynamic disorder and a scattering from phonons, 
we can increase the effect of the temperature to the DOS 
using the eq. (dzl)- 

Furthermore we shall consider the effect of the real part 
of the self energy, from eq.(25) we have: 


A 

riicr^ 


a(A_ - A+) -/3r2 In 


^2 1 


A_A. 




Al + Al 


(29) 

2 -1 

and using the previous result In ^ and 

being multiplied by F^ thanks to eq.(27) it will be ex¬ 
ponentially small and it can be neglected. Defining 
A_ — A I = ^5 We obtain: 

^ cx. 


A = 


A^ 


Al 


(30) 


which is an amazing result, because it introduces a lin¬ 
ear dependence of the Dirac point with the variance of 
the disorder and it moves to negative energies increas¬ 
ing We deduce that the disorder can shift the Dirac 
point. Taking into account the low disordered nature of 
Graphene, the effect it will have a physical meaning just 
for a value of the variance up to cr^ ^ 0.1. In particular 
also in this case if we consider just the static disorder 
we can shift and fix the position of the DP, fixing the 
variance of the impurities in the system and considering 
again eq.([T7| we can also move the DP linearly with the 
temperature in a continuous way, managing the position 
of the gap of the system with the temperature. In partic¬ 
ular if we substitute the numerical values of the constants 


in eq.(30) we obtain: 


And in the particular case of t = 0.2t we obtain: 


(31) 


(32) 


We can see those effects of eq.(30) and eq.(27) in fig. 
where we observe both the shift and the non zero value 
of the DOS. 

In the previous section we have seen that the effects 
induced by the asymmetry are quite relevant from an in¬ 
tellectual point of view, but when we use the real param¬ 
eters of graphene we notice a small effect. We have also 
seen that the effects the disorder caused the shift of the 
DP and they increase the effects of the asymmetry. For 
this reason we can perform a comparison between the nu¬ 
merical results of the tight binding integration with and 
without interaction. In particular we can analyze the ef¬ 
fect with ^ — 0 and at high temperature with a minimum 
number of electrons in a band equal to no = 2 and the 
maximum equal to n = 4. 

We can first analyze the fi = 0 case. In the DP the disor¬ 
der increases the effect of the thermal doping, as we can 


Dirac point shifting 



FIG. 2. Continuous shifting with increasing disorder. It is a 
numerical calculation with in the range from 0 to 0.1, with 
an increment of — 0.02. The picture is in hopping units 
with t = 0.2t 


see from fig.([^. In particular we can obtain an analyti¬ 
cal approximation of the amplified thermal doping fitting 
the red curve in the figure, obtaining the function: 

f(T) = 2 + aT^!'^ (33) 


where a = (1.532 ± 0.005) in hopping units. We can see 
a comparison between eq. and eq.p3| also in fig.®, 
with a net increase. In particular the expression eq.( |33[ ) 
may be helpful for a quantitative comparison with the 
experimental datas. Furthermore if we consider the ther¬ 
mal doping at room temperature (T = 300A') we obtain 
the following result with and without disorder: 


half filling(T = 300A:): 

3 g2 109 varucles 


• half filling(T = 300Ar) with disorder cr^ = 0.09: 
n - 2.27 10i° 


So with the disorder we can increase the thermal dop¬ 
ing by one order of magnitude. It is interesting also to 
analyze the effect of the disorder at high temperature. 
In particular in this case we can analyze the number of 
particles as a function of temperature and also the chem¬ 
ical potential and as we can see in fig.Q in both cases 
the disorder can increase the effect of the thermal dop¬ 
ing. These properties might be important not only for 
the physics of graphene, but also for interesting applica¬ 
tion as a semiconductor. In particular we did not report 
the low temperature limit because it is less important 
for application. Furthermore considering that the single 
particle solution is it valid only for a restricted value of 
temperatures, by eq. ([TtI), we know the low temperatures 
correspond also to low values of the disorder, so for hav¬ 
ing significative differences between the disorder and the 
ordered case we should go beyond physical values of the 
cr2 ~ O.I. 
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H=0 numeric 
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Dirac point 


V. COMPRESSIBILITY 


The compressibility of the Fermi gas is a well known 
property which can be measured by means of a Scanning 
Tunnelling Microscopy (STM). In fact thanks to the STM 
we can have maps of two dimensional conductors and 
semiconductors and we can also measure the DOS. We 
can define the compressibility as: 


K = 



(34) 



FIG. 3. Termal Doping in the Dirac point with and with¬ 
out disorder calcnlated nnmerically (top) and comparison be¬ 
tween the analytic solutions without disorder and a fit with 
fixed disorder (bottom). The picture is in hopping units with 
t = 0.2t and no = 2 


T=0.06 



Far from the Dirac point 



T/t 

FIG. 4. Numerical results for the doping with Ghemical po¬ 
tential with fixed temperature (top) and with temperature 
with a fixed chemical potential (bottom). The picture is in 
hopping units with t — 0.2t,no = 2 and /ro = 0.6f 


where there is a very high sensibility for measuring the 
charge fraction, as reported irP^. If we know the volt¬ 
age between the tip and the sample we can relate it to 
the inverse compressibility as At4 = as reported 

iiP3! 

for graphene. In particular in this way it is possible 
to analyze the effects of the disorder on the system. Fur¬ 
thermore it might be easier to measure the effect of the 
impurities and the e-ph interaction on the compressibil¬ 
ity than the thermal doping and for a Fermionic system 
we can use the Einstein relation between the conductiv¬ 
ity and the compressibility, i.e. a = eD ( j j where D 


is the diffusion coefficient linked to the disorder of the 
system, as reported irP^. 

We can use the functional form of eq.(27) for calculat¬ 
ing the compressibility as K = , following the Seitz 

theorem, as reported in Maharl^. In particular inside 
the SCBA approximation we have seen that the DOS of 
the system is a local function independent of k and in 
the case of the static disorder it is also indipendent of 
/i, so deriving by /i and dividing the DOS in two parts 
respectively for electrons and holes, we shall obtain: 


r n/oc / ^9fiu; + p) 

+ ‘‘K hi 

Jo dp 


(35) 


that is an exact expression for a disordered system. In 
particular if we consider just the static disorder, we will 
have to include just the first order correction and we 
can neglect the Hartree contribution, which shifts the 
energy to a density dependent value Eij(n), dependent 
on p only, as reported in eq.(15). In particular the prop¬ 
agator of eq.(I2) will present the additional correction of 
E//(n). The total compressibility of the system can be 
written as: 


X = (l-^) (36) 

where p = p — po — T,h, Po = is the DP shift. 

So we can write = —oy and the compressibility will 
be: 


Xo 

I - axo 


( 37 ) 
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where x is the total compressibility and xo is the com¬ 
pressibility with just the Hartree correction, we recognize 
in eq.(37) a structure similar to the RPA. In particular in 
the next subsection we shall analyze the simple impurity 
scattering correction and then the e-ph contribution. 


A. Hartree contribution 


and solving the integrals we will obtain linear corrections 
in temperature and cubic corrections in the chemical po¬ 
tential as: 

={Pr + h + 

® , ,,3 («) 

+ (/3B-/30(2ryic,g(2) + 5^) 


We can calculate the compressibility from the impu¬ 
rity scattering considering the expression of the DOS re¬ 
ported in fig.([^ and fitting the expression with an easier 
form of the DOS, expressed as a polynomial expansion 
and with a non zero value in the DP, as: 


DOS[uj) = ± (38) 

in this way the compressibility will be made by a con¬ 
tribution from the first order and another one from the 
second order, considering that the derivative by p, is equal 
to a derivative by w with an opposite sign and with LjR 
distinction between electrons and holes. In particular 
also in this case there will be a behavior near the DP 
and another one far from the DP. In particular for the 
low temperature result we can use again the Sommerfeld 
expansion for the first order as: 


/ TT^ „ , 

d^ni ~ DOS{y. ) -f -^^OS ifi )T^ -- 

~ 

^ V ' 6 (r2 + (a2^')2) 


and for the contribution of the second order we will have: 


d^>n 2 = 


± PR/L 


+ 2 xT^i) (4( 


and so the total low temperature compressibility is: 

- + 


, - 2 2 p2 

d,n= + —T"-^ 

^ V ^ 6 (r2 + (a,/,')2)i 

±/3fl/L + 


We notice that for T —)■ 0 the compressibility reduces to 
the DOS of the system. The high temperature limit will 
be quite similar, but with a further correction from the 
integrals of the second order, because near the DP we 
shall consider both the contributions from electrons and 
holes: 


= 

= Pl j dx{ - ^^^){{xTY + {^i)‘^ + 2xTh)-{- ^^ 2 ) 

+ Pr J + + 2 xTn ) 


regarding the second order, we can use a box integration 
in the interval [—1,1] dividing the integrand by 2. In this 
way the total compressibility near the DP will be: 


x = i^ = ^l(^ + T)\l 1 + + ry- 


4T 




+ -m+ 


+ 


r2 


4aT L 
+ (Pr + Pl){ 


asinh(^(/r -|-T) - asinh(^(/r -T)) 
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+ {PR-PL){2T^P log (2) + ^) 


(44) 


B. Fock contribution 


We have seen that the interactions in the system can be 
modeled like a RPA correction. In particular we can see 
that calculating a we will be able to increase the effect of 
the interaction in the system. In particular the constant 
a, using the values of the deformation potential of the 
previous section, will be: 

a = 2|ffq^oP-D(q = 0,w„ = 0) 
a ~ 1.3(t units) 


and so using the eq. (371 we will be able to increase the ef¬ 
fect of the asymmetry and disorder to the compressibility. 
We can see the results in fig.(VB). Furthermore consid¬ 
ering the temperature dependence of the compressibility 
at the Dirac point, we can speculate on the temperature 
dependence of the conductivity with a = assum¬ 


ing a weak dependence of the diffusion coefficient from 
the temperature. In particular we obtain the following 
fitting function: 


/(T) = AT (46) 

where we notice a linear dependence and A = (10.95 ± 
0.07). We can observe the numerical result with the fit¬ 
ting function in figure |VB[ where we notice that the 
points near the zero, which are the most important one, 
because near the values used in the experiments, they do 
not present an exact linear dependence, but it is not easy 
to verify this dependence just numerically, especially near 
the zero. 
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FIG. 5. Quantum capacitance as a function of the chemi¬ 
cal potential at low temperature (top) and high temperature 
(bottom) with a comparison between numerical and analyti¬ 
cal calculation. In both pictures kf, = 1, the energy scale is in 
hopping units t = 0.2t. 


T(K) 

16 32 48 64 



T/t 

FIG. 6. Quantum capacitance in the DP for low temperature 
computed numerically for t = 0.2t in hopping units and in 
Kelvin. We notice the linear dependence, but with a deviation 
for a temperature under l&K. 


VI. CONCLUSIONS 

We can finally analyze the results in the paper. In par¬ 
ticular we have seen we took inspiration for our work from 
an experimental paper where they measures an anomal 
conductivity and furthermore introduced a new experi¬ 


mental set-upP. We tryied to make some basis for those 
unexpected phenomena. In particular we have seen that 
the effect of the e-h asymmetry at low energy induced 
an asymmetry in the DOS, which is perfectly symmet¬ 
ric without t . Moreover we have seen that the number 
of particles of the system are functions of the tempera¬ 
ture just throughout /3, which is the e-h correction. This 
kind of correction allows us to move the Fermi level of 
the system with temperature, but even though the prop¬ 
erty is interesting for the physics of graphene, there is a 
small effect on real material. This aspect has induced the 
study of disorder and its combination with the e-h asym¬ 
metry. The impurity scattering with the SCBA approach 
allows us to increase the thermal doping at room temper¬ 
ature by one order of magnitude and considering also the 
correction to the DOS, at hxed energy, we found strong 
effects in the DP. In particular thanks to the disorder we 
can shift the DP increasing the variance of the disorder 
and with the dynamic disorder and the e-ph interaction 
this shift will be temperature dependent. This property 
is especially evident in the calculations of the interacting 
DOS and the compressibility. 

We have seen that exists an experimental measure of t . 
It may be possible to measure these properties with the 
Scanning Tunneling Spectroscopy (STS) and the Angle 
Resolved Photoemission Spectroscopy (ARPES). In par¬ 
ticular the STS can measure the DOS of the system, 
for measuring the DOS and its asymmetry and with the 
ARPES is possible to measure the band structure of the 
material and in particular the e-ph coupling in the sys¬ 
tem. There are some previous studies about the effects 
of irradiations on graphene sample^^, using an electron 
beam with 3QKev, registering a shift of the DP of ^ 20V 
for graphene on Si 02 and this is probably caused by the 
strong interaction with the substrate otherwise they also 
saw DP shift for suspended graphene, up to ^ 0.16U 
and it suggests that disorder may induce this property. 
This is a possible way for verifying these effects and also 
manipulating the material adhoc. In fact the STM can 
manipulate single charges and in this way we may be able 
to choose and shift the insulating point according to our 
necessities. Furthermore as possibles perspectives we no¬ 
tice that we have not considered yet the electron-electron 
interaction, and very recently they studied the Coulomb 
interaction in association to the disordeil^, and also the 
combined effect with the spin for considering short and 
long (RKKY) magnetic interaction and t he q uantum spin 
hall effect, because the Kane-Mele modeP^l is similar to 
the t correction. Moreover the last year an experimen¬ 
tal group from StanforcP^ discovered a three dimensional 
Dirac material as the NasBi, with a band structure pre¬ 
senting Dirac cones. New research prospectives may be 
to study the e-h asymmetry for this kind of materials, 
originated by a next nearest neighbors in a tight-binding 
model. Probably thanks to the new geometries it may 
possible to find new thermal and magnetic effects, or 
strong impurities interactions, arriving to new techno¬ 
logical applications in electronics. Furthermore really 
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recently emerged the possibility of making honey-comb 
lattices of cold atoms^, where they obtained a DP band 
structure and a quantized conductivity, so it may be pos¬ 
sible to apply similar calculations to those system. 
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VII. APPENDIX 1-DOS 

Using eq.Q the expression for calculating the DOS is: 

p{E) = J dq5{E - E±{q,e)) (47) 

Now using the Dirac delta function properties: 


r*+oo 


N 


dx5{f{x)) = ^ 


2=1 


' —OO 


l/'(u)| 


(48) 


we directly obtain the q integral with the delta argu¬ 
ments: 


9i.2(e,6») = 


±vf ± ^v)-AA{9y 


and so: 


P(e) = 






dO 


2A{e) 

9i(e,^') 


^Vf + 2A{6)qi[t,9) 


(49) 


(50) 


Considering our interest for the DP we can expand our 
solutions up to the second order of e and using directly 
eq.Q where with the q{e, 9): 


gi(e,6i) 




m 




(51) 


and so thanks to eq.(49) and eq.(50) we can expand the 
integral: 


P(e) = 




p27r 


(27r) 


d9 E \ 2> 


A{9) 


(52) 




and solving the eq.(52), the final result will be equal to 
the volume of the unitary cell for the integrand function 
multiplied by 27r, because there is just one angular de¬ 
pendence and it is zero. Finally the expression of the 
DOS without a degeneracy factor g = 4 is; 


P(e) = 


3-\/3a^ |e| 

47r v'j 


± 


27v^a^ / t 


Stt 


T -3 (53) 


like it has been reported in the text. 


VIII. APPENDIX 2-DOPING 

A. High temperature limit 

Inside this regime we can use the result of eq.([^, but 
with p, y 0 and \p\ T, obtaining some corrections 
to the previous behavior. Taylor Expanding the Fermi- 
Dirac function up to the second order: 

fie, p, T) - fie, 0, T) + fie, 0, T) ~ fie, 0, T)+ 

( 54 ) 


-f 


df 

dp 


p -I- 


d^f 

dp^ 


and using it inside eq.([^ we will have eq.([^ plus other 
correction with a p dependence: 


nip,T) = 

r~p(e) f/(e,0,r) 


and so: 

nip,T) - u(0,r) = 


p/1 


\d^f] 

dp 

/i -h 

fi —0 

_dp^ _ 




\dn 


\d^n 

L ( 

dp 

{E—O 

dp^ _ 





(56) 


we can write the derivatives in p as derivatives in e as 

0 ^ consequent dependence of 
the functions from the energy. Considering the different 
contributions relative to electrons and holes and the first 
and second order: 

nip,T) - n(0,T) = 




\df] 


\d^n 

1 

.de. 

fi—0 

de^ _ 


- fi—O 


\df] 


\d^n 

.de_ 

11 + 

/x=0 

de‘^_ 


/2 = 0 



= niip,T) E n 2 ip,T) 


(57) 


where niip, T) considers just contributions from first or¬ 
der and n 2 ip,T) from second order. Using again the 
substitution ^ = x and with a change of the sign of the 
integrals regarding p-(e): 


niip,T) = p {2aTh E 2l3T^l2) 


where every is given by: 


T — T 

^71 — 


n 




2 

and in our case U = log 2 I 2 = ^. 
In a similar way we obtain: 


(58) 


(59) 


n 2 ip,T) - n(0,r) = 

9/(e)^+~ 


P±(e)- 


de 


p deiaE2ye)i-^) 


(60) 


and so neglecting exponential small contributions: 

n 2 ip,T) = -p^aIoE 2 pTh) (61) 

and putting together everything we obtain what has been 
reported below. 
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IX. APPENDIX-3 


X. APPENDIX-SPECTRAL PROPERTIES 


In particular we can find an expression for the self 
energy considering some properties of the involved in¬ 
tegrals: 


/„(z,A)= / de- 


(62) 


This is the full expression of the complex logarithm 
with the corrections: 


where e = uj(h = 1) and uj — E(e) = z. They present 
some recursive relations such as: 


n In z - e 


/o 

n—1 

de - 


cA 


de- 


and so: 


/„(z, A) = A) - 


A" 


with the properties: 


-^+) — ) 


(63) 

(64) 

(65) 


( 66 ) 

(67) 


so using eq. ([^ and the eq. (65) we arrive to an expression 
for the self energy of the system: 


S(£) 

71 * 0-2 


= aGiz) + PF{z) 


where 


G{z) = A_ — A+ -I- z In 


F[z) = zG{z) - 

In particular using z = e we have: 

^ = (a + M 

A_ — A_|_ -f 07 In 


(z-f A_)(z-A+) 
A^ -f A?. 


( 68 ) 

(69) 

(70) 


(w -I- A_)(a; — A+) 


- (71) 


logz^ = 


r(cc) 


log ((w — A(w))^ — r(a;)^) -I- 2i arctan — , sgn(a;) = 
^ ^ w — A(w) 

log Pi -I- 2i9i 

— log z -I- A_ = 

1 *^7r 

- - log (07 - A(o;) + A_)2 - r(a;)2 - 7Sgn(o;)(- + 


arctan 


r(o7) 


07 — A(o7) + A_ 


) = 


- 2 log ^2 - id2 

- log z - A+ = 

- t log (07 - A(o;) - A+)2 - r(07)2 - *sgn(07)(-| + 


arctan ■ 


r(o7) 


' — A(o7) — A+' 


- 2 loS P3 - idz 


(72) 


and r(o7) and A(o7) are the the expressions reported in 
eq.pl. 
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